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We examine the quantum tunneling process in Bose con- 
densates of two interacting species trapped in a double well 
configuration. We discover the condition under which parti- 
cles of different species can tunnel as pairs through the poten- 
tial barrier between two wells in opposition directions. This 
novel form of tunneling is due to the interspecies interaction 
that eliminates the self- trapping effect. The correlated mo- 
tion of tunneling atoms leads to the generation of quantum 
entanglement between two macroscopically coherent systems. 

PACS numbers: 03.75.Fi, 03.65.Ud, 74.50.+r 



Quantum tunneling of macroscopically coherent sys- 
tems is an intriguing phenomena well known in the con- 
text of Josephson junction effects in superconducting 
electronic systems. For super-fluids consisting of neu- 
tral particles, detailed investigations of tunneling is aided 
by the recent realization Bose-Einstein condensation of 
atomic vapor in a well controllable environment. In- 
deed, recent experiments have successfully demonstrated 
quantum tunneling for condensates confined in an ar- 
ray of optical potentials [0,1). One prominent feature of 
tunneling in Bose condensates is the nonlinear dynamics 
arising from the interaction between atoms. Quite re- 
markably, for single-component condensates trapped in 
double- well configurations, previous studies have indi- 
cated that a self-trapping mechanism can suppress the 
tunneling rate significantly by increasing the atom-atom 
interaction strength ||||||||]. 

An interesting extension of the tunneling problem in- 
volves Bose condensates of two interacting species (Fig. 
1). The main issue is how the interspecies interaction af- 
fects the tunneling process, and particularly the quantum 
coherence as the two condensates mix together. Previ- 
ous studies of the general properties of two-component 
Bose condensates have emphasized the important role of 
the interspecies interaction, which leads to novel features, 
such as the components separation (0,|), cancellation of 
the mean field energy shift 0], and the suppression of 
quantum phase diffusion fl(i|j . However, the investigation 
of the influence of interspecies interaction on tunneling 
dynamics has only just begun [pd| , ^2| . 

In this paper we present a novel tunneling mechanism 
for a two-component condensate trapped in a double- 
well (see Fig. 1). The atoms of the component A(B) are 
initially prepared in the left (right) potential well. We 
discover the condition under which the interspecies in- 
teraction can eliminate the self-trapping effect and thus 



enhances the tunneling significantly. Such an enhanced 
tunneling originates from the correlated quantized mo- 
tion of the two condensates. We also show that atoms of 
different species tunnel through the barrier as correlated 
pairs in opposition directions, i.e., a form of quantum en- 
tangled tunneling. Therefore tunneling serves as a mech- 
anism to build up a strong correlation among atoms of 
different species, and this leads to the generation of quan- 
tum entanglement between two multi-particle systems. 

The configuration of our double- well system is sketched 
in Fig. 1. Our focus in this paper is the quantum dynam- 
ics beyond the mean field description. An exact many- 
body description is difficult even for single-component 
condensate problems. The usual method to capture the 
essential physics is based on the two-mode approxima- 
tion in which the evolution is confined by the left and 
right localized mode functions associated with the respec- 
tive potential wells 0,0]E |Jlj§ ■ Such an approximation is 
valid when each potential well is sufficiently deep so that 
higher modes of the wells essentially do not participate 
in the dynamics. 




FIG. 1. A sketch of the double- well tunneling system. Ini- 
tially, the atoms of the component A are trapped in the left 
side and the atoms of the component B are trapped in the 
right side. 

In the two-mode approximation, the system is modeled 
by the Hamiltonian (U = 1), 
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Here the subscripts L and R respectively denote the lo- 
calized modes in the left and right potential wells. Since 
there are two modes available for each component, the 
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model in fact consists of four mode operators. We use 
aj and fet (j = L, R) to denote the creation operators of 
the component A and B respectively. The parameters f2, 
K a (nb) and k describe the tunneling rate, self- interaction 
strength of the component A{B) and the interspecies in- 
teraction strength respectively. 

To gain insight of the quantum correlation develop- 
ing in the tunneling process, we first consider the ex- 
actly solvable case with only one A atom in the left 
well and one B atom in the right well. In this case the 
system is spanned by four basis vectors: |1, Q)a|1> 0)_b, 
|1,0)a|0,1} b , |0,1}a|1,0} b and |0, 1) A |0, 1) B , where 
\p, q)s denotes the state with p atoms of species S (S = 
A, B) in the left well and q atoms of species S in the right 
well. The eigenvalues and eigenvectors of H can be found 
straightforwardly. In the regime where the interspecies 
interaction is sufficiently strong such that k 3> O, the 
state vector evolves as 

|\fir(t)) = e-^^+^/^^cost^ll, 0} A |0, 1) b 

+ismu} Q t\0,l) A \l,0) B ]+O(n/n). (2) 

In writing Eq. (^|) we have defined luq = Q 2 /2k as an 
effective tunneling frequency. Because of the strong in- 
teraction between the atoms, the probability of finding 
both particles in the same well at any time t is neg- 
ligible (of order f2 2 /« 2 ). The tunneling motion of the 
two atoms are anti-correlated in the sense that the atom 
A and the atom B always move in opposite directions. 
Such an anti-correlated tunneling motion gives rise to 
quantum entanglement between the two atoms. At time 
t = (n + l/4)7r/w0) ( n = integers), the state is a form 
of Bell's state that is maximally entangled in the two- 
particle two-mode subspace. 

Now we examine the multiple atoms case. In or- 
der to facilitate the discussion, we assume the number 
of particles are the same for the two components, i.e., 
N a = Nb = N, and the condensates shares the same 
interaction strength i.e., n a = Kh = k- The latter con- 
dition is a good approximation to 87 Rb condensate of 
atoms in hyperfine spins states \F — 2, my = 1} and 
\F = 1, mt = — I), which share similar scattering lengths 
||. However, we emphasize that these assumptions are 
not crucial, we shall relax these conditions later in the 
paper. 

We shall limit our study to the An^Ntt regime where 
the nonlinear interaction is dominant. As before we con- 
sider the initial condition in which all atoms in the com- 
ponent A{B) are localized in the left (right) potential 
well. The general form of the state vector at time t is 



y/(m + l)(N-m) Cn,m+1 



%N t • 



E^=o E™=0 Cn,m{t)\n,N 



given by: |^(t)) = 

n) A \m, N — m) B . The amplitudes c„ )TO (i) are governed 
by the Schrodinger equation according to the Hamilto- 
nian dll): 
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yJ{n+l){N 



+y/m(N - m + l)c n 
+k [(n + m- N) 2 ] C n ^ n 



(3) 



with the initial condition c„ !TO (0) = <5 n ,JV<5m,o- 

Let us first present the typical results obtained from 
the numerical solutions of Eq. (3). In Fig. 2, we show 
the particle number difference W = (a L aL — Qfl a fl) of 
species A between the two wells as a function of time. 
The occurrence of tunneling is revealed by the decrease 
of W. At longer times W approaches zero, therefore the 
numbers of A atoms in the two potential wells are roughly 
equalized. We emphasize that the nonzero interspecies 
interaction is responsible for the tunneling to occur. If 
the two species do not interact with each other (i.e., k = 
0), then a sufficiently strong self-interaction Kj > Nil 
(j = a, b) can suppress the tunneling almost completely 
by the self-trapping effect . 




FIG. 2. Particle number difference of the component A be- 
tween the two potential wells as a function of dimensionless 
time (Q. 2 / ' n)t with N = 5 and k = n a — Kb = lOfi. The 
dashed line is the solution based on the effective Hamiltonian 
in Eq. (6). The inset shows the overlap probability P of 
finding the system belonging to the degenerate set of states 



+ ^Jn{N - n + rjcn-i. 



To understand the mechanism of the two-component 
tunneling, we note that the state vector (3) is a superposi- 
tion of (A^+l) 2 states of the form \n, N—n) A \m, N—m) B . 
However, only a few number of them are actively in- 
volved. This is intuitively understood because \n,N — 
n) A \m, N — m) B have different energies for different val- 
ues of n and m, and only those with energies closer to 
that of the initial state are more accessible. 

The difference in energies is significant in the regime 
4k 3> Ntt considered here. We find that the states 

\4> 9 ) = \q,N-q) A \N-q,q) B (q = 0, 1, 2, N) (4) 

are approximately degenerate energy eigenvectors of H 
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(to zero order in fi) , and the energies of states other than 
\4> q ) are higher than that of \<p q ) by an amount of order of 
k or higher. Since the initial state |\&(0)) is |^o), is 
mainly a superposition of \<f> q ) at any time t according to 
the energy argument above. We find that this is indeed 
the case. To provide a numerical evidence of our finding, 
we show in the inset of Fig. 2 the overlap probability 

N 

defined by P(t) = \ (cj) q \ &(t)) | 2 . For the parameters 

used in this figure, the set of \4> q ) contributes more than 
90% of Our further numerical tests suggest that 

P —> 1 in the limit NQ/k -> 0. 

We remark that all \<fi q ) have the same number of par- 
ticles (component A plus component B) in the left po- 
tential well, and the same also holds for the right well. 
Therefore P(t) w 1 implies small fluctuations of total 
particle number in the each potential well. In the case of 
Fig. 2, we find that the fluctuation of total particle num- 

1 /2 

ber in the left well (ANl) is much smaller than (Nl) , 
i.e., a sub-Poissonian distribution. 

The time-dependent problem now is significantly sim- 
plified because the evolution of the condensates mainly 
involves the set of degenerate states \<p q ). For those 
states that do not have the same energy as \<fi q ), they 
act as intermediate states that are rarely populated. We 
may eliminate such intermediate states with the adia- 
batic following method which is well known in quan- 
tum optics. The last line in Eq. (3) indicates that the 
states with n + m = N + r (r =integer) has the diag- 
onal term r 2 K that can be interpreted as the energy of 
the states if the tunneling interaction Q were switched 
off. Since k is a large parameter here, the states with 
n + m = N ± 1 have a much different energy than that 
of the states \4> q ) with n + rn = N. Any transition (due 
to Q) from the manifold n + rn = N to the manifolds 
n + m = N ± 1 must quickly return to n + m = N. In 
other words, the states with n + m = N ± 1 are interme- 
diate states that are hardly occupied at any time. The 
amplitudes associated with the n + m = N ± 1 manifold 
are: c„ + i^-n, Cn-i,N-n> Cn,N-n+i and c n ,N-n-i, and 
they can be found approximately by adiabatic following 
under the condition 4k^>N£1. With these approximate 
amplitudes, the Schrodinger equation of c„ N-n is given 
by: 

n 2 

ic n ,N-n~--^-(n + l)(N - n)Cn + x N-n-1 

ZK 

n 2 

- — n(N -n+ l)Cn-\ ,N-n+l 
ZK 

n 2 

[n(N-n) + N]cn !N - n . (5) 

K 

This corresponds to the Schrodinger equation governed 
by the effective Hamiltonian: 

fl 2 ( 

H cS = -— (a^aR&i&Jj + a L a^ L b R + a\a L a\a R 

+b[b L b R b B )j , (6) 



apart from a constant term proportional to the total 
number of particles. The effective Hamiltonian H c f[ is 
an approximation that captures the essential tunneling 
mechanism in the 4k ^> NO, limit. We have tested the 
validity of the effective Hamiltonian numerically For ex- 
ample, in Fig. 2 the dashed line, which is obtained from 
the evolution based on H Q g, agrees well with the exact 
numerical solution. 

The physical picture becomes transparent according to 
H c g. The interaction term a^b^aLbfi suggests that every 
time when an atom A moves from the left well to the right 
well there must be an atom B moves from the right well 
to the left well. This explains why the small fluctuations 
in the total particle number in each potential well. The 
reverse process is described by a^b^aB.bL in H c g. In 
other words, the atoms A and B must be moving in pair 
in opposite directions during the tunneling process. 

It is worth noting that H c g can be cast in the form 

H eS = -^-(K + K_+K_K + ) (7) 

where K + = a) L aR + b^ L b R and K- = a^ajj + bi,b R 
satisfy the angular momentum commutation relation: 
[K+,K_] = 2K Z , [K X ,K±] = ±K±, where K z = (a{a L + 
b^jpi, — aJjOR — b R bu)/2. Therefore K± and K z are anal- 
ogous to collective spin operators, and the eigenvectors 
and eigenvalues of H c ff can be solved analytically using 
angular momentum algebra. This shares a common fea- 
ture in spinor condensates where the spin mixing problem 
can be solved in a similar fashion Jli| . The nonlinear in- 
teraction between collective spins is the key to generate 
nonclassical correlations between different spin compo- 
nents, and particularly, the notion of multi-particle en- 
tanglement has been discussed in the context of squeezed 
spin states |l^Jl^,|l^| . We show here that quantum entan- 
glement between two multi-particle subsystems (^4 and 
B) can be achieved in the double-well tunneling process. 

To this end we quantify the degree of entanglement 
between the two species based on entanglement entropy: 
E = -tx{pA\ripA) = -tr(p B lnp s ), where p A and p B 
are reduced density matrices of the respective subsys- 
tems, i.e., p A = tr bPab and p B = tr aPab with p A B = 
\^(t)) {9(t)\ being the density matrix of the whole sys- 
tem. A disentangled state (for example the initial state 
above) has zero entanglement entropy. The more entan- 
gled the systems, the larger the value of E is. As an 
illustration, we show in Fig. 3 how the entanglement is 
established in time for various particle numbers. As the 
time increases the value of E increases until a saturated 
value is reached. Since there are N + 1 degenerate states 
\4> q ) mainly involved in the evolution, E w In (N + 1) is a 
maximum if all \(j> q ) have equal contributions to the state 
of the system. In the case of Fig. 3, the value of E can 
reach as high as 90% of the value ln(iV +1). 

Finally we would like to address the conditions for 
the entangled tunneling to occur. Our discussion above 
has been restricted to the simplest symmetric situation: 
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N a = Nt — N and K a = Kb — k, in order to illustrate the 
essential mechanism under the condition 4k ^> Nil. The 
same analysis can be performed to examine the general 
situations. We have examined the system with unequal 
particle numbers (Nb — N a ) = D > 1 and unequal cou- 
pling strengths n a = k + 6, Kb = k — S with \S\ <C k. We 
find that if the tunneling strength Q is sufficiently weak 
or the self interaction strength k is sufficiently strong 
such that 



A\k(D - 1) - [N a (D - 2)|J|]/2| » QN b . 



(8) 



then the system mainly evolves among the states \n, N a — 
n) A \N a - n, N b - N a + n) B , where n = 0, 1,2, ...N a . In 
other words, the tunneling under condition (8) is char- 
acterized by entanglement generation term a} R b\aLbR as 
before. However, we point out that the tunneling are gen- 
erally less efficient for the cases with non-zero 5 and D. 
This is because the states |n, N a — n) A \ N a — n,Nb — N a + 
n) B are not as degenerate as that in the symmetric case 
with 5 and D are both zero. We also remark that the con- 
dition (8) does not apply to the special case Nb — N a = ±1 
in which some of the states in the n + m = N, and 
n + m = N ± 1 manifolds are accidentally degenerated. 




(Q, z /K)t 



FIG. 3. The time dependence of the entanglement entropy 
E for the different particle numbers: iV = 2, N — 10 and 

iV = 50 With K — K a — Kb — 100. 

To conclude, we have presented a novel mechanism of 
double-well tunneling involving Bose condensates of two 
interacting components, based on the two-mode approxi- 
mation model in the strong coupling regime. We find that 
the interplay of intraspecies and interspecies interactions 
permits a set of energy degenerate states, a small tun- 
neling coupling can push the system to 'explore' through 
these degenerate states and thus result in a substantial 
tunneling not limited by the self trapping effect. The 
most interesting feature is the strongly correlated tun- 
neling motion. We have derived an exactly solvable ef- 
fective Hamiltonian to capture the mixing dynamics. In 
addition, we have shown that high degree of quantum 



entanglement between two macroscopically coherent sys- 
tems can be achieved. 
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